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Abstract. We study topological properties of quasi-particle states in the non-Abelian quantum 
Hall states. We apply a skein-theoretic method to the Read-Rezayi state whose effective theory is 
the SU{2)k Chern-Simons theory. As a generalization of the Pfaffian {K = 2) and the Fibonacci 
(K = 3) anyon states, we compute the braiding matrices of quasi-particle states with arbitrary spins. 
Furthermore we propose a method to compute the entanglement entropy skein-theoretically. We 
find that the entanglement entropy has a nontrivial contribution called the topological entanglement 
entropy which depends on the quantum dimension of non-Abehan quasi-particle intertwining two 
subsystems. 



1. Introduction 



Since the discovery of the fractional Hall effect and the Laughlin wave function, it has been 
extensively studied the fractional statistics, i.e., anyon (see e.g. Ref. |59|). The effective theory 
of the Laughlin wave function for u = 2^+1 the U{1) Chern-Simons (CS) theory (SOi 16411 . 



and those quasi-particles as excited state are understood as the Abelian anyon. Contrary to these 
fractional quantum Hall states, the u = 5/2 Hall state experimentally observed in Ref. l60l remains 
mysterious (see e.g. Ref. l5l|) . One of theoretical candidates for this state is the Moore-Read 
Pfaffian state n44\ . which is related to the BCS ground state wave function of the p + ip supercon- 
ductor such as Sr2Ru04 11241 15211 . The effective theory of the Pfaffian state was identified with the 
SU (2)2 CS theory 11211 . and the non-Abelian property and the dimension of the Hilbert space of 
the many-quasi-particle state have been studied from this viewpoint [,461 . 



As a generalization of the Pfaffian state as the non-Abelian quantum Hall state, proposed is the 
. This state with i^' = 3 is expected to describe the u = 12/5 quantum Hall 



62L and also pointed out is a relationship with the rotating Bose-Einstein 



Read-Rezayi state [|53 
states reported in Ref. 

condensates [8]. Originally the Read-Rezayi state was constructed by use of the parafermion 
theory 16311 which is closely related to the SU {2)k CS theory. Algebraic structure of the Read- 
Rezayi state was analyzed by use of the quantum group Uq{SU(2)), and the non-Abelian property 
as the braiding matrices of quasi-particle states was studied nSw . 
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The CS field theory plays an important role in the quantum topological invariant of links and 3- 
manifolds. Witten constructed the quantum invariant of 3-manifolds from the Chern-Simons par- 
tition function with gauge group G based on conformal block ||45^, and he further identified quan- 
tum polynomial invariants such as the Jones polynomial [29] and the HOMFLY polynomial [|l7|1 
as the expectation value of the Wilson loop of the CS theory [16 IN . His construction in terms of 
path integral is not mathematically rigorous, and the invariant of 3-manifold was later constructed 
combinatorially by Reshetikhin and Turaev [|54ll . and it is called the Witten-Reshetikhin-Turaev 
(WRT) invariant. 



Recently from the viewpoint of the quantum computation (see e.g. Refs. [33, 1471) . the non- 
Abelian quantum Hall states have received renewed interests. As a fault-tolerant quantum com- 
puting system, Kitaev proposed to construct quantum registers topologically i34l l36ll. Though his 
toric code uses a spin system on 2-dimensional honeycomb lattice as a non-Abelian quasi-particle 
system (see also Ref. |40 where topological phase of trivalent graph, sometimes called string net 
or spin network, is discussed), the quantum Hall system is another candidate for non-Abelian 
statistics (see Refs. Isl @, 20 for recent reviews). In such topological quantum registers, unitary 
operations are constructed from braiding of non-Abelian anyons. It was shown [|l3L 11411 that the 
braiding in the SU {2)k Chem-Simons theory can efficiently approximate, i.e., quantum compile, 
any quantum computing gate; n>3-braid group for > 3, K 7^ 4, 8 and n>5-braid group for 
i^" = 8 are capable of universal quantum computation. 

One of our purposes in this article is to study the braiding property of non-Abelian quasi- 
particles in the SU{2)k CS theory. By use of explicit form of correlation function in Ref. |38 . 
Ardonne and Schoutens studied braiding matrices Jll for 4-quasi-particle states. It becomes dif- 
ficult to derive braiding matrices for many-quasi-particle states by their method. We rather use 
the Temperley-Lieb algebra or the skein theory, which is elementary and combinatorial, to derive 
the braiding matrices. Although, a naive application of the skein theory does not give a unitary 
representation. We show explicitly how to construct the unitary braiding matrices by use of skein 
theory. See also Ref. 32, where proposed is a different method to construct the unitary operators 
from the skein theory by modifying a weight of trivalent vertex. Our method is essentially equiv- 
alent to that of Ref. Qo, once we know both a construction of the colored Jones polynomial based 



on representation theory of the quantum group and a quantum group structure of the conformal 
field theory (see e.g. Refs. |23i |30|). So most of methods in this article might be well known, 



except our new result is a computation of the entanglement entropy (see Refs. 1 for recent re- 
view) in the framework of the skein theory. In a topological phase such as non-Abelian states, 
it was discussed pSi . 1391] that the entanglement entropy has a topological part, which is called 
the topological entanglement entropy, and that it is expected to detect a topological order [|58n . 
Some numerical and analytical studies have been done for the quantum dimer model in a trian- 
gular lattice [18], for the fractional quantum Hall states of the fermionic Laughlin type [0], and 
for the quantum eight- vertex model [|48|]. Here we show that the entanglement entropy of the state 
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1^) = J2jPj l^i)A ® with V. |p,f = 1 is given by 



'3 



which includes the topological entanglement entropy 

j 

Here dj denotes the quantum dimension of quasi-particle which connects subspaces A and B. 

A content of this article is constructed as follows. In Section [21 we review the Read-Rezayi 
state. We shall discuss a relationship with the Zk parafermion CFT, the SU{2)k Chern-Simons 
theory, and the SU{2) topological quantum invariants. We further explain basic facts about the 
Temperley-Lieb skein theory following Ref . 31. In Section [3] we outline a method to compute 
the entanglement entropy of quasi-particle states skein-theoretically. We define the topological 
entanglement entropy from the von Neumann entropy. In Section |4] we construct a unitary rep- 
resentation of braid operators on spin-1/2 state in the SU{2)k CS theory. The SU{2)2 theory 
corresponds to the Moore-Read Pfaffian state. We shall further compute the entanglement en- 
tropy of these states, and we show that it has a topological entropy which depends on the quantum 
dimension of quasi-particle. Section [5] is devoted to studies of spin-1 quasi-particle states in the 
SU (2) K theory. A case of = 3 is known as the Fibonacci anyon model. We construct the 
bases of quasi-particle states, and give explicit forms of the braiding matrices. We also calcu- 
late the topological entanglement entropy. In Section |6] we study a correlation function which is 
composed from different species of quasi-particles as another candidate of topological quantum 
registers. We give an elementary treatment of skein theory up to Section|6]as best as we can. Later 
in Section |7] we shall give a general setup for many-quasi-particle states in terms of the quantum 
6j symbol. The braiding matrices and the topological entanglement entropy are discussed by use 
of the quantum 6j symbol. The last section is devoted to concluding remarks and discussions. 

2. Non-Abelian Quantum Hall State, CFT, and Quantum Invariants 
2.1 Parafermion Theory and the Read-Rezayi State 



The Read-Rezayi state is constructed by use of the parafermion theory [|2lL 16311 (see also 
e.g. Ref. 10). The Zk parafermion theory is equivalent to SU{2)k/U{1) theory, and we use the 
primary fields $f„ with —£< m < i and i E {0,1,..., K} whose conformal dimension is 

^ _ i{i + 2) 
™~4(ir + 2) AK 

Here m denotes the U(l) charge defined modulo 2 K, and parameter i parametrizes a spin-£/2. 
The primary fields Gm with m G Z/2 of the SU {2)k theory is 

G^^^ = $Lexp(A^^ (2.1) 
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- m 

= ^L+2K = ^m^K- The fusion rule follows from the SU {2)k theory as 



where (p is the free massless Bose field. Further for the field $^ we have i = m mod 2 and 



min{e+e',2K-e-e') 

*ix<,= <+rn' (2.2) 

e'=\t-e\ 

In terms of the primary fields the vacuum sector / is ${] = <l>^ = /, and the parafermions 
the spin fields a^, the dual spin fields and neutral fields si are respectively defined by 



(2.3) 



ae = ^[ 
H = ^-i 
ee = K 

and we mean ipl = ipK-e, and aj = ax-i- Note that h^^ = = h^^ = 2k^k+2) ' ^^'^ 

K+2 • 

It is noted that the field for K = 2 becomes the Majorana fermion tp satisfying 
{ip{z) ^(w)) = (z — w)^^, and the Pfaffian state is identified with the correlation function 



(ipizi) ■ ■ ■iI){zn)) ll44ll . Correspondingly the wave function of the n-quasi-particle states can be 



constructed from the correlation function (cri(wi) ■ ■ ■ ai{wn) ipizi) ■ ■ ■iI){zn)). This function was 



studied in detail ll46ll . and the dimension of the wave function for 2 n-quasi-particle states are 
identified to be 2""^ 



The Read-Rezayi state 115 311 is a generalization of the Pfaffian state, and the wave function is 



defined as correlation function of the primary fields in the parafermion theory; 

^(Wl, ...,Wn]Zi,...,ZN) = ((Ti(wi) • ■ ■(yi{Wn)iJl{zi) ■ ■ ■ i)l{zN)) 

N n 

X n (-^^ - '^^•)'''^^ n n - ^j-)* n - ^.o'^™^ (2.4) 

i<j i=l j=l i<j 

Here M is an odd or even integer depending on whether the particles are fermions or bosons, and 
the filling factor becomes 

K 

V 



KM + 2 

The fractional power in the second line of (I2.4|) was added to cancel the singular power which 
comes from the operator product expansion of so that the wave function is non-singular. We 
should note that the Read-Rezayi state (12.41) without quasi-particles {n = 0) is an exact ground 
state of the Hamiltonian with (K + l)-body interactions; 

n = V ^ 6'^{zi^ - z,^)5'^{zi^ - Zi^) ■ ■ ■6'^{zij^ - Zij^^J (2.5) 

il<i2---<iK+i 
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So to Study the action of braiding of quasi-particles on the wave function (12.41) . we need the 
braiding actions on the correlation functions of parafermion fields. By use of correspondence (12.11) 
with the SU{2)k theory, we see that the braiding matrices coincide with that of the SU{2)k 
theory. 



2.2 Quantum Invariant from the Chern-Simons Theory 



In Ref. |6li the physical interpretation of the quantum knot invariant such as the Jones polyno- 
mial and the HOMFLY polynomial [17] was given. The Chern-Simons action with gauge 
group G on 3-manifold M. is defined by 



M 



(2.6) 



where Khz. coupling constant K E Z, and A is the gauge connection on Ai. Witten introduced 
the topological invariant of M as the partition function 



Z{M) 



VAe 



(2.7) 



This invariant was later mathematically rigorously defined by Reshetikhin and Turaev [1541], and it 
is called the Witten-Reshetikhin-Turaev (WRT) invariant. Simpler construction was given later 
in Ref. |4l| based on the Conway-Kauffman skein relation (see also Refs. yjj, l42r) . 



Not only the invariant of the 3-manifold Ai, the Chern-Simons action (|2.6I) gives the quantum 
invariant for link Cm Ai. We define the Wilson loop expectation value by 

1 



Ri,...,Rj 



VAW'^^^{A)...W^^-{A)e 



(2.8) 



Z{M) 

where we assume the link C has components /Ci , . . . /C^, and denotes the Wilson loop operator 
along a component /C with a representation R of gauge group G; 

W^{A) = Ttr V exp ^ a (2.9) 

K. 

where P is a path ordering. When the gauge group G = SU (2), this correlation function becomes 
the colored Jones polynomial at the root of unity. 



2.3 Combinatorial Construction of SU (2) Quantum Invariants 



Hereafter in this article we only study a case of the gauge group G 
as a S'-transformation of the conformal block, we use [|22|] 

(a + l) (6 + 1) vr 



ab 



sm 



SU{2). In SU{2) case, 



(2.10) 



K+2 V K+2 
In actual computation of the quantum invariants such as Z(Ai) and Wr{]C), useful and the most 
elementary method, at least to the author, is to use the skein theory. Here we briefly review a 
skein-theoretical construction of the SU (2) quantum invariant (see, e.g., Ref. yjj, |42|) . 



6 



K. HIKAMI 



The Jones polynomial is characterized by the skein relation; 



A 



) ( -H^-' 



(2.11) 

DVJ Q =dD (2.12) 

Here A is a parameter related to the quantum group deformation parameter q = A^, and 

d = -A^-A"'^ (2.13) 

Projection of knots and links is denoted by D. 

In constructing the colored Jones polynomial of links and the SU(2) WRT invariant of 3- 
manifolds, fundamental tool is the Jones-Wenzl idempotent. Throughout this article, we employ 
a standard notation in diagrams D\ an integer n beside an arc means that there exist n copies of 
that arc. Like (|2.1 II) an arc without integer denotes a single arc, i.e., = 1 is often omitted. Then 
the Jones-Wenzl idempotent is depicted as a blank square, and it may also be recursively defined 
as 

n + 1 



A 



n-l 



A, 



n - 1 



Here A„ is the colored Jones polynomial for unknot; 

An = 

A2n+2 A-2n-2 

= f-ir " 




A^ - A-2 

We have Aq = 1, Ai = ci, and others are recursively defined by 

^n+l = d An — A„,_i 

Note that we have 

Ax + y + Z + l Az-l = Ar^^z Ay^z " Ay Aj 

One notices that the Jones-Wenzl idempotent satisfies foUowings; 



m 



n + 2 









m 




n 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



where n > m. 
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In the SU{2)k theory, we have A„ 7^ for \ < n < K, and ^k+i = 0; namely the parameter 
A satisfies 

^4(X+2) _ ^ (2.19) 

This condition means that integer n attached to each arc should he 1 < n < K. Here we 
choose il5|| 

vri 



A = i exp 



2{K + 2) 

from physical requirement which will be explained below. Under this condition we have 



sin tt) 



and especially we have 



d = Ai = 2 cos 



TT 



K + 2 



We introduce a trivalent vertex by 




where 



a + b — c 



m 



c + a — b 



b + c — a 



n 



2 2 
Lidices in trivalent vertex (12.231) should fulfill an admissible condition; 

a + 6 + c = mod 2 

a + b > c, 6 + c>a, c + a >b 

a + b + c<2K 

Skein relation (12.1 II) gives 



-A' 




-A- 



Generally we have a twist formula; 




A 



ab 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



where 



(2.27) 
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2.4 Trivalent Diagram as Wilson Line of Quasi-Particles 

The diagram introduced above can be interpreted physically as follows. We regard a single arc 
as a Wilson line of quasi-particle with spin- 1/2. The Jones-Wenzl idempotent which has index n 
means the projection to spin-r2/2 space, and the arc labelled by n corresponds to the Wilson line 
of quasi-particle with spin-n/2; 



□ : spin-n/2 quasi-particles 



Then the trivalent vertex (12.231) denotes the usual fusion rule; quasi-particles, (pa and (pf,, with spin- 
a/2 and spin-6/2 fuse to spin-c/2 quasi-particle 0^- One sees that the admissible condition (12.241 ) 



of the SU {2)k CS theory has a correspondence with the fusion channel (|2.2I) Dol . In general the 
fusion rule is written as 

0aX05 = J]iV,^fe0c (2.28) 



Verlinde formula [|570 shows that the the fusion multiplicity can be written in terms of the S-matrix 
as 

a 

We recall that the ^-matrix for SU (2) is given in (12.101) . The quantum dimension da of quasi- 
particle with spin a/2 is the largest eigenvalue of the matrix Na, and in the SU{2)k case we 
have 

da = ^ = Aa (2.30) 

•J 00 

The non-vanishing correlation function of quasi-particles should fuse into the vacuum section. 
For example, the correlation function of two quasi-particles with spin- 1/2 will vanish unless two 
fuse into the vacuum. Diagrammatically we denote this 2-quasi-particle state by 

|qp)= W (2.31) 

Two endpoints of arc indicate existence of quasi-particle at spatially different points, and a con- 
necting arc denotes that they fuse to the vacuum sector. Naturally a dual state is defined upside 
down as 

(q~p|= (2.32) 
We can compute norm of the state by connecting two endpoints of (qp| and |qp), and (12.121) gives 

(qp|qp)= O ^^-^^^ 
As we have set a parameter A as (|2.20l) . we have (i > as a level K is K > 1. Then we have the 
normalized physical state, which is schematically written as 

|qp) = i= W (2.34) 
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Entanglement is one of distinguishing properties of quantum mechanics, and it receives much 
interests in recent studies of quantum information science. One of measures of entanglement is 
the entropy (see, e.g. , Refs. IH for recent review). Namely when we assume that the state is a pure 
bipartite state |\1'), and that the system is divided into two sub-systems A and B, the Schmidt 
decomposition assures that we can write 

\'^)=J2p^ \^3)a®\<I^j)b (3-1) 
j 

where ^ and ^ are orthonormal states in subspace A and B, 

^(V^J^fe)^ = ^j,k B'\^Mk)B = ^hk (3.2) 



and ^ \pj\^ = 1. Then the entanglement entropy is defined by 



■'J 



Sa = -^IPjf'^oglpjf (3.3) 
j 

We note that 

Sa = Sb (3.4) 

Alternatively the entanglement entropy is defined as the von Neumann entropy 

Sa = - TiA (pa log Pa) (3.5) 
Here pA is the reduced density matrix (Alice's density matrix) 

PA = TTBP = TTBm m) (3.6) 

where the density matrix p is 

P=|^)(*| (3.7) 
In fact when we use the replica trick 16. 261 

Sa = - lim ^ {TtaP^} (3.8) 

n-*i on 

we obtain (|3.3I) . 



In the non-Abelian states, it is expected that there exists an effect from a topological order 135 . 



3911 . Those papers discuss that such effect appears from a dependence on a length of boundary 
between A and B. As our model is purely topological and pure gauge theory, we rather define the 
topological entanglement entropy by 

ST° = - Tr A {pA log Pa) + ^ \Pj I ^ log \pj \ ^ (3.9) 

j 

We shall compute the topological entanglement entropy (13.91) of the non-Abelian quasi-particle 
states in the following. 



10 



K. HIKAMI 



We explain how to compute the topological entanglement based on the skein theory. In our 
formulation, the state is schematically written as 



(3.10) 



where each end of arcs denotes quasi-particles at spatially different points. There exists a sum of 
trivalent graphs insides an ellipse, which is a Wilson line describing how quasi-particles fuse to 
the vacuum sector. Correspondingly the density matrix (13.71) is depicted as 



(3.11) 



We now divide quasi-particles into two groups A and B. We mean that the owner of each quasi- 
particle is Alice or Bob. As the Alice's reduced density matrix pA (13.61) is given by taking traces 
in Bob's space, it is depicted as follows by connecting ends of Bob's arcs; 



PA 




(3.12) 



where in the last expression, A and B denote quasi-particles belonging to Alice and Bob respec- 
tively, not spins of quasi-particles. Applying (13.81 ). the von Neumann entropy (13.51) is computed 
by connecting n copies of the reduced density matrix p^; 




(3.13) 



Above trivalent graphs are computable by the skein theory, and we can obtain explicitly the topo- 
logical entanglement entropy from (13.91) . 



4. Quasi-particles WITH Spin 1/2 

We study the braiding operations on the correlation function of many-quasi-particles with spin- 
1/2. The Bratteli diagram of spin- 1/2 quasi-particles is depicted in Fig. \T\ One sees that the 
number of quasi-particles must be even so that the correlation function is non- vanishing. In the 
case of 2-quasi-particle state, we have only a unique path which starts from (0, 0) and ends at 
(2, 0) in the Bratteli diagram in Fig.[Tl which corresponds to the state (12.341) . Generally we have 
several paths from (0, 0) to (2 n, 0). The number of paths denotes that of the fusion channels, and 
it is the dimension of the Hilbert space of 2 n-quasi-particle states. We see that, for the SU (2)2 



theory, it is 2" ^, which denotes the dimension studied in Ref. |46 
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We shall study the action of the braid operators on the 2 n-quasi-particle states. The non- 
Abelian property of quasi-particles is that we have the nontrivial representation of the braid oper- 
ators. Hereafter (jj means the braid operator acting on the i-th and i + 1-th particles as 



i 1 + 1 



(4.1) 



Here i and i + 1 mean the i-th and i + 1-th positions from the left respectively, and these should 
not be confused with spins of quasi-particles. The operators Ui satisfy the Artin braid relation 





(4.2) 



Easy is to see that 



for |i — j I > 2. We also use the twist Of, 



UiUj — (Tj Ui 



(4.3) 



4.1 4-Quasi-Particle States 

We study the braiding effects in the 4-quasi-particle states in the SU {2)k>2 theory. A dimension 
of the Hilbert space of 4 quasi-particles is two, and we may choose unnormalized bases as 

|0)= U W 

ID- ^ 
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Correspondingly dual bases (0| and (1| are defined upside down. The norms of these diagrams are 
computed as 



(0|0)= =d' 





where we have used 



(4.4) 



We see that they have positive norms, d? > Q and A2 > 0, under a parametrization (|2.20l) . It is 
noted that, when K = 1, we have A2 = 0, and that 



1 ) becomes a null state due to that it does 



not satisfy the admissible condition (12.241) . We also see that 



(0|1) 



by making use of the following identity, which is a special case of (12.181) : 







Thus we indeed have two normalized physical states 

1 



10)^^ U U 



|1) 



which are orthonormal; for i.j G {0, 1} we have 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



We study actions of braid operators on states {|0), |1)}. Here we use the F-matrix defined by 
We see that 

/ i 1 \ 

(4.9) 



■pii 



1 
d 

1 - 



1 

_L _1 

which can be computed by use of a recursion relation of the Jones-Wenzl idempotent (|2.14|) . 
Actions of ai and can be computed easily by use of (|2.25l) . For instance, we have 



V 



-A \j \j 
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To get an action of a2, the F-matrix (I4.8h will be useful. We have 



13 



0-2 



/ 1 y I \ 



ill 

11 



pill 
*11 



ill 

11 




v 




-A^ 

A, 

-A^ 

A, 



Actions of other operators can be computed in the same way. In summary, the representation of 
the braid operators on the space spanned by {|0), |1)} (14.61) are as follows; 



p(ai) = p((T3) 



-A- 



A 



d Ad 



(4.10) 



y/dP^ 1 
Ad A^d 



which are unitary under (12.201) . We can check directly that these representations satisfy the Artin 
braid relation (14. 2L and that 

P ((Ti (T2 Ci^ 0-2 0"1 ) = A^^ ■ 1 



Computation of the representation of the twist (14.31) is straightforward. As all quasi-particles 
are spin- 1/2, the representation of 9i is same for all i, and we have 



Pie) = -A' 



(4.11) 



It is remarked that these representations are computed in Ref. |2 
functions as the 4-point correlation functions of the WZW model 



from an explicit form of wave 

3. 



4.2 6-Quasi-Particle States 



We next study 6-quasi-particle states. We restrict to the case K > 2. The SU (2)2 case will be 
discussed separately in the later section. As can be seen from the Bratteli diagram in Fig. [U the 
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dimension of the Hilbert spaces is 5. We set 

1 



|00) 
|01) 
|10) 

111) 



d \/d 



w w ^ 



1 



{<p - 1) 

d {d^ — 1) 



(4.12) 



- 1) (rf2 _ 2) y 2 V 2 y 

We have defined these states based on Ref. '20'. See that (12.201) proves positivity of norm of 
each diagram, and by the same computation with a case of 4-quasi-particle, we see that these are 
orthonormal bases 



where i,j E {00, 01, 10, 11, C}. 



We study the action of the braid operators {(Xi, . . . , 0-5}. Actions of cti, (T3, and 0-5 are easily 
computed by use of (12.261) . Actions of a2 on |00) and |10), and of 0-4 on |00) and |01), can be 
given in the same method with 4-quasi-particles. For the action of others, we may use 





(4.13) 



Here the F-matrix is given by 



^12 
12 



d2-l 



1 



(4.14) 



which is derived by the recursion relation of the Jones-Wenzl idempotent (|2.14|) . Then, for exam- 
ple, we compute as follows; 



0-2 KJ 




(F!l)„(-^-') +(F!l)u-4 



ill 

11 



+ F 





A 

A-'^ 
A 
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In the last expression, we further apply (14.131) to get 

^^^^X^ =mr)n y ^.y y.y.y 

Other actions can be computed in the similar methods. 

To summarize, the representation of the braid operators on the 5-dimensional space 
{|00),|01),|10), 111), are given by 

p(cTi) = diag {-A-\ -A~^, A, A, A) 



( -'- 









d 



Ad 





Ad 



V 



Wd 



1 

Ad 
Ad 










1 

Ad 



d 
Ad 



\ 



Ad 



Ad 



p((T3) = diag (-A '\A,A,-A 3, A) 



' Ad 



d 
Ad 




\ 



1 



















A3 


1 


d 


Ad 


1 




Ad 


d 






Ad 


Ad 






Ad 



Ad 



(4.15) 



p(a5) = diag (-A \A,-A \A,A) 
One can check that these satisfy (14.21) , and 

P (cTi 0-2 0-3 0-4 0-5 0-5 0-4 0-3 0-20-1) = A~^ ■ 1 
Note that the representation of the twist (14.31) is same with before, and we have (14.1 II) . 



4.3 Topological Entanglement Entropy of Spin- 1/2 Quasi-Particle States 

%^ 2- Quasi-Particle State 

As the first simple example of explicit computation of entanglement entropy, we study the 2- 
quasi-particle state (|2.34l) . We have the state |\1') = |qp), and the density matrix p is depicted 
as 

1 KJ 

p = - (4.16) 
d ^ 

We assume that both Alice and Bob have one quasi-particle. Alice (resp. Bob) has the left (resp. 
right) quasi-particle, and we denote it as A = {1} and B = {2}. As a trace Tr^ means connecting 
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arc ends of Bob's quasi-particles (13.121) . we have 

P{i} = Tr|2} p = i 



(4.17) 



Then from (13.131) we obtain the entanglement entropy as 



and the topological entanglement entropy (13.91) is identified with 

5;°7 = %} = logc? (4.19) 

This result is due to that Alice and Bob are intertwined by a Wilson line of spin- 1/2 quasi-particle 
whose quantum dimension is d (I2.30|) . 

%^ 4-Quasi-Particle States 

We next consider the 4-quasi-particle states (|4.6I) . which spans two-dimensional space. We take 
a state as 

m=Po \0)+Pi |1) (4.20) 

where |j5o| + = 1. We first assume that the first and the second quasi-particles from the 
left belong to Alice, and that the remaining third and fourth quasi-particles to Bob; we denote 
A = {1, 2} and B = {3, 4}. Connecting ends of Wilson lines of Bob's quasi-particles as in (|3.12l) 
and using (14.41 ) and (14.51) . Alice's reduced density matrix is schematically written as 

2 1 W 2 1 V 

P{i,2| = bo| d ^+\P^\ ^ (4.21) 

It is straightforward to compute (13.131) . and we get the von Neumann entropy as 




boriogbor - bil'log I 1 (4.22) 



where we have used (14.41) and (14.51) . To extract a topological entropy, we notice that the state (14.201) 
has the orthogonal bases 2} l'^i){3 4} (13.21) thanks to (14.41) and (14.51) . As a result, the 

topological entanglement entropy (|3.9I) is identified with 

s';;^;^ = \p^flog{d'-l) (4.23) 

This result is interpreted as follows. The state |0) is not topologically entangled between Alice 
and Bob, or they are entangled by the vacuum sector with quantum dimension do = 1. Although, 
in the state 1 1) Alice's quasi-particles are intertwined with Bob's through the Wilson line of quasi- 
particle with spin-1 whose quantum dimension (12.301) is = A2 = (i^ — 1. 
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We make similar computations for the 4-quasi-particle states (14.61) when Alice owns the second 
and the third particles while Bob has remainder; A = {2, 3} and B = {1,4}. Following (|3.12l) 
the reduced density matrix is given by 

|2 



P{2,3} 



rf2 



\Pl 



+ 



PoPI+P*oPi 



+ 



\P1 



d?-l 



(p-l 



□ + 



d ^ 



(4.24) 



where 



Po 



Po , 



d 



Pi = Po 



1 



(4.25) 



Pi 
d 



Here, in the second equality, we have substituted a definition of a trivalent vertex (12.231) and 
have used the Jones- Wenzl idempotent (I2.14|) . Note that |po|^ + = 1- Using (|2.18l) the 
entanglement entropy (13.131) is computed as 

5'{2,3} 





- bii^iog 



rf2 - 1 



- bo I ^ log IpoI 



(4.26) 



To extract the topological entanglement entropy, we should note that bases in (14.61) do not give 
orthogonal bases, |V'j){2 3} I0i){i4}' Alice's and Bob's spaces (13.21) . To remedy it we 
rewrite the state (|4.20l) by use of the F-matrix (|4.8I) as 



1^) 



P'd 



(4.27) 



which has a form of (13.21) . i.e., the states \i'j)^2 3} \'Pj){i4} orthonormal due to (|4.4I) 



"J/{1,4} 



and (14.51 ) (see, e.g., (17.81 )). Thus the topological entanglement entropy (13.51 ) is given by 



ritopo 
{2,3} 



Ipii'log (d^-l) 



(4.28) 



When we look at the expression (14.271) , we see that Alice's and Bob's quasi-particles are not 
topologically entangled with a probability \pof while they are intertwined with a probability \pif 
by the Wilson line of quasi-particle with spin-1 whose quantum dimension is d2 
the expression (14.271) is comparable with (14.201 ) replacing pa with pa. 



^2 



1. Thus 



We leave it for readers to check that 



ritopo 
•^{1} 



QtOpO 

{2,3,4} 



hgd (4.29) 

which supports our interpretation because Alice and Bob are intertwined by spin- 1/2 quasi- 
particle (consult computations in Section |74l ). In view of (14.231) . (14.281 ). and (14.291 ). the topological 
entanglement entropy depends on what quasi-particles Alice and Bob have in the state |^). Even 
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if the state is same, there exists a possibility that a different quantum dimension appears depending 
on a quasi-particle which intertwines Alice and Bob when we change owners of quasi-particles. 



4.4 Pfafflan State 

The Moore-Read Pfafflan state is described by the SU (2)2 CS theory, and we have A = i e^^^^ 
and d = We then have a representation of the twist Qi (14.31) as 



(4.30) 



The braid matrices (|4.10l) for 4-quasi-particle states (|4.6I) become 

p((T2 



e 8' 



1 i 



V2 V 1 



(4.31) 



For 6-quasi-particle, the state \C) amongst (14.121 ) is not admissible due to (12.241) ; we can see 
that norm of the diagram for |C) reduces to zero, and it is unphysical, null state. Then from (|4.15l) 
the braiding matrices on space spanned by bases {|00), |01), |10), |11)} in (14.121 ) are given by 



e"'/*^ dia^ 


li- 


-1,- 


-1 






(l 





i 


o\ 


g-|7ri 





1 





i 









1 









i 





V 


e"'/® dia^ 






i, 


-1) 




/I 


i 





o\ 


g-|7ri 


i 


1 
















1 


i 




\^ 





i 


V 


e"'/^ dia^ 











(4.32) 



These representations are studied in Ref. [l9|. Also see Ref. |28| where studied are the braiding 
matrices of vortices in a p-wave superconductor. 



Concerning the topological entanglement entropy, we see that the quantum dimensions (|2.30l) 
are given by di = \f2 and ^2 = 1, and that only spin- 1/2 quasi-particles contribute to the topo- 
logical entanglement entropy. Our results, (|4.23l) and (|4.28l) . indicate that the topological en- 
tanglement entropy vanishes in the 4-quasi-particle states (14.201) when both Alice and Bob have 
two quasi-particles. Note that it does not vanish when Alice has one or three quasi-particles 
(see (14:291) ). 
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5. QUASI-PARTICLES WITH SPIN-1 

We study the braid group operations on the correlation functions of quasi-particles with spin-1. 
The fusion rule of spin-1 gives the Bratteli diagram as in Fig. |2l As in the previous section, the 
number of paths from (0, 0) to {n, 0) is the dimension of the conformal block of ra-quasi-particles. 
See that, different from the spin- 1/2 case, we have the non-vanishing correlation function even 
when n is odd. Both the dimension of 2- and 3-quasi-particle states is one, and the normalized 
2-quasi-particle and 3-quasi-particle states are respectively given by 

In Iqpg) two endpoints of arc labelled 2 indicates that there are two quasi-particles with spin-1 in 
spatially different points. Correspondingly in Iqpg) three ends of arcs denote that there are three 
quasi-particles in spatially different points. 



6- ► ► ► 




2- ► ► ► ► ► 




1 i 5 J s f 



Figure 2: Bratteli diagram for spin-1 for SU{2)7 theory 



5.1 4-Quasi-Particles States 



We first assume K > 3. The SU (2)3 theory will be discussed later. A dimension of the Hilbert 
space of 4-quasi-particles with spin-1 is three as can be read from the Bratteli diagram in Fig. [2l 
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1 



rf2 - 1 2 2 



|1) 



d 



2 2 2 2 

1 ^ 



(d2 - 2) 



(5.3) 



2 2 2 2 



Vd^ - 3 rf2 

A parameter (12.201) supports positivity of norm of each diagram, and the recursion relation (|2.14l) 
of the Jones-Wenzl idempotent proves the orthonormality; 



The F-matrix which we use in this section is 

2 2 



H 



p22 
22 



2 2 

X 

2 2 
2 2 



(5.4) 



where 



-p22 
■^22 



1 

(P-2 
d{cP~l) 



d 

1^ 



1 - 



d2-2 



(5.5) 



\ 1 _ d."^ __d(d^-3£+l) 1 I 

(d2-l)2 (d2_l)(rf2_2)2 (d2„l) (d2_2)/ 



This identity follows from 



2" ' "2 
j2 2 



1 
d 



3 2 



A2A3 



which can be derived after some algebra. 



Actions of the braid operators on bases {|0), |1), |2)} (15.31) are computed similarly. Represen- 
tations of cTi and can be given from (12.261 ) straightforwardly, and that of 0-2 can be given once 
we use the F-matrix (15.41) . Explicitly they are written as 

p(ai) = p(a3) = diag {A-^-A-\A') 



( 



P(0^2) 



d2-l 

\ A4(d2_i) 



1 1 



Vd4_3^2_|_x 

A4 (d2-l) 



1 



d2-2 



AS (d2_2) 



d^-3d2+l 
d2-l 



(5.6) 



' d4-3d2+l 1 , 

i (rf2-2) Y ((i2_l) (d2_2) / 



See that the braid relation (14.21) is fulfilled, and that 

p(ai (72 0-3 0-3 0-2 O"! 
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The twist (|4.3I) is follows from (12.261) because all the quasi-particles are spin-1; 



21 

(5.7) 



5.2 5-Quasi-Particles States 



In SU {2)k theory with K > 5, the dimension of the Hilbert space of 5-quasi particles is 6. We 
choose bases as 



|01) 
1 10) 

111) 
|12) 

|21) 

|22) 



d 



2 2 2 2 2 



LV 



A2 VA3 

fl 22222 

2 2 2 2 2 



A2 

d' 



22222 



d 



(5.8) 



v/A^ 



4 2 4 

2 2 2 2 2 



VA3A4 4 2 

22222 

VLV 



dA. 



A. A, 



^4 ^5 4 4 

Here |a6) means that intermediate states have spins-a, b as can be seen from (|5.8I) . and we give 
generalization later by use of the quantum 6j-symbols. We have normalized these states to be 
= 6ij, where dual states {i\ are defined using upside down trivalent graph in (15.81) with the 
same normalization factors. 

To obtain the representation of the braid operators, we need another F-matrix beside (15.41) : 



H 

2 4 

H 



\ 



Tp22 
*24 



\ 2 4 / 



2 4 
2 2 

X 



(5.9) 



where 
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1 - 



A3 ^ 

(f J_ 
(As)^ A3, 



(5.10) 



A twist formula (12.261) gives representations on spaces spanned by 

{|01),|10),|11),|12),|21),|22)}; 



p(ai) = diag {A-\ -A~\ -A~\ -A~\ A\ A') 
p((T4) = diag {-A-\ -A-\ A\ -A-\ A^) 



(5.11) 
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Using above mentioned F-matrices, we obtain 



/ At. 



















yl4 





As 
Aa A3 










18 _d_ 

^ A3 



/SA7 



751 





A^3 



Aa A3 





P((^3) 








V 

We can see that these satisfies (14.21) . and 

P ((Ti 0-2 0-3 0-4 0-4 0-3 0-20-1)= A^'^^ 

The twist is represented as (15.71) . 



\ 







A3 





d , 


















AS 


A* 


VAI 








Aa 




A-iAa 


A4 


As 


d /A4 










Aa A3 


AS A3 Y Aa 




rf /a7 


d 








Aa 


AS A3 Y Aa 


AI6 Aa A3 











A^^ 
A3 


VdK^o 

A3 











VdA; 
A3 


AS A3 



(5.12) 



5.3 Topological Entanglement Entropy of Spin-1 Quasi-Particle States 
%^ 2-Quasi-Particle State 



We start from the 2-quasi-particle state Iqpa) (15.11) . Here we suppose that the first quasi-particle 
belongs to Alice and the second to Bob; A = {1} and B = {2}. By use of (|3.12l) the reduced 
density matrix is depicted as 

2 



1 



from which (|3.13l) gives 



= = log {<f - 1) 



(5.13) 



(5.14) 



Recalling the topological entanglement entropy (14.191) of the two spin-1/2 quasi-particle states, 
we see that d is replaced with — \ which is nothing but the quantum dimension of the spin-1- 
quasi-particle who intertwines Alice and Bob. 
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%^ 3-Quasi-Particle State 

In the 3-quasi-particle state |qp3) (15.21) . we suppose that the owner of the first and the second 
quasi-particles is Alice while the third is Bob; A = {1,2} and B = {3}. Taking a trace, i.e., 
connecting both ends of the Wilson line of the third particle (13.121) . we get 

2 2 
2 2 

By definition (I3.13|) . we get the topological entanglement entropy 

■^E} = ^{1,2} = log [d^ - 1) (5.16) 

This agrees with the logarithm of the quantum dimension of the spin-1 quasi-particle which inter- 
twines Alice's and Bob's quasi-particles. 



%^ 4-Quasi-Particle State 

We compute the entanglement entropy for the 4-quasi-particle states defined by 

|^)=Po |0)+pi \l)+P2 |2) (5.17) 

Here bases \a) are defined in (15.31) . and \pq\^ + \pi\^ + \p2\^ = 1. As the first example, we group 
four quasi-particles as A = {1,2} and B = {3, 4}. Through (13.121) we get the Alice's reduced 
density matrix as 

2 2 2 2 

The entanglement entropy (13.131) can be computed as 

%2} = -\Poflog\po\'-\p,f\og (^^j - \p,f\0g (^ ^.Xi^^i ) ^^-^^^ 

As bases (|5.3I) have a form of (13.21) due to (|7.8I) . the topological entanglement entropy (13.91) is 
given by 

^Im} = bli'log - 1) + b^l'log (t/^ - 3t/2 + 1) (5.20) 

We see the appearance of quantum dimension d2 = d'^ — 1 and ^4 = — 3 rf^ + 1 of quasi-particles 
which respectively intertwines Alice and Bob in the state |1) and |2). 
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When we assume that A = {2, 3} and B = {1, 4}, it is convenient to change bases using the 
F-matrix (15.41) . We have another expression of (15.171) : 

^22 ^ 2222 



where 



1 1 Vd^-3d^ + l 



1 rf2_3 1 /(;4_3^2^1 



Vci^ - 3 ^2 + 1 1 /(;4_3^2 + i 1 

P2 = Po 15 z Pi -j: \ 15 h P2 



rf2 - 1 rf2 _ 2 V _ 1 ' (^^2 _ 1) (^2 _ 2) 

See that |po|^ + + |p2|^ = 1- By the same computation, we obtain the (topological) entangle- 
ment entropy as follows; 

5{2,3| = - IPOI^ log IPol^ - \Plf log (^^) - log (^ a^Xi^^l 1 ^^-^^^ 

^lls} = l^^i I ' log - 1) + |P2 1 ' log (rf^ - 3 + 1) (5.24) 

which is consistent with our interpretation of the topological entanglement entropy; when we look 
at (15.211) . we find that parameters pi and pz respectively denote probabilities that the spin-1 and 
spin-2 quasi-particle intertwines Alice and Bob. 



5.4 Fibonacci Anyon 



The spin-1 quasi-particle in the SU{2)3 theory is known as the Fibonacci anyon [|49|1 . This is 
because the dimension of the Hilbert space of n-quasi-particle with spin-1, i.e., the number of 
paths from (0, 0) to (n, 0) in Fig. |2l coincides with the Fibonacci number. For instance, among 
three states of 4-quasi-particles (15.31) . the state |2) is not admissible (|2.24l) ; the norm of the diagram 
for |2) vanishes, and it is a null state. Setting A = ie'^'/^^ and d to be the golden-ratio d = 
2 cos(7r/5) = ^-y^ = r in (15.61 ) we obtain 2-dimensional representation of the braid operators on 
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{|0),|l)}as 



= p(ct3) 



oh 




(5.25) 



Pid) 



65' 



which are given in Ref. |49| assuming the unitarity condition. 



A dimension of the Hilbert space of 5-quasi-particle states becomes 3; the norms of diagrams 
in the states |12), |21), and |22) vanish due to A4 = 0, and these states are not admissible (12.241) . 
Note that only |22) is forbidden for a case of i^T = 4. Then on bases {|01), |10), 1 11)}, the braiding 
matrices (|5.11|) and (|5.12l) reduce to 



diag (^e-5^',e5^',e 




\ 

/el-' 





65" 







C5 1 



\ 

eh 



(5.26) 



p(ff4) = diag(eS'",e-S",et") 



On the topological entanglement entropy (13. 9L we recall that the quantum dimension of quasi- 
particles with spin- 1/2, 1, and 3/2 are respectively given by di = r, c/2 = r, and ^3 = 1. Thus 
the topological entanglement entropy always takes a form of |p^|^logr where \Pt\^, satisfying 
< \Pt\^ < 1, is a probability that Alice and Bob are intertwined through quasi-particle with 
spin- 1/2 or spin-1. 



6. QuBiT FROM Different Species of Quasi-particles 

In previous sections, we have only considered states which are constructed from quasi-particles 
with same spin; spin-1/2 and spin-1. In this section we shall propose another candidate of topo- 
logical qubit which is composed from 4-quasi-particles with different spins. Namely we set 
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1 < n < K — 1, and study states constituted from two spin- 1/2 quasi-particles and two spin- 
n/2 quasi-particles. We set bases as follows; 



Oi) 
O2) 
O3) 
O4) 

O5) 
Oe) 



— In In 

n - 1 

- « 1 In 

1 y_y 



1 « In 



n - 1 

1/7 « 1 



An 



n - 1 
« 1 n 1 



A„ 
1 



1 



n - 1 



111) 
II2) 
II3) 
II4) 

II5) 
lie) 



7^ 

V ^n+l « + 1 



^ n I I n 

V n + 1 



In n I 



7^ 



/J 1 n 1 



V^n+l n + l 



An A„_|_i 

A„ A„+i 



11 « H 



M n 11 



(6.1) 



Dual states {ia\ are upside down of \ia) with the same normalization factor. A condition (12.211) 
proves that these states are physical states, i.e., norm of each diagram is positive. Normalizations 
are chosen so that 



{ia\3a) = 



fori, j G {0, 1} and 1 < a < 6. 

To compute explicit representation of braid operators, we prepare the F-matrices. Due to the 
definition of the idempotent, we find 



" \ /I 

1 / \n 



(-1) 



n + l 



1 - 



1 



In + 1 
I ^"^n 



dA„ 



17, - 1 

Ti + 1 



(6.2) 



(6.3) 
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Using these F-matrices, we can compute the actions of braid operators on 
{|0i), |li), IO2), II2), ■ ■ ■ lOe), lie)} are computed. The actions of ai and 0-3 follow from (12.261) as 



P(^i) 



/ 


Ri 











\ 


Ri 


























Ri 














Ri 























R2 





\ 














T-1 ; 

R3/ 







Ri 








°^ 











Ri 








Ri 




















Ri 


























R3 





\o 














R2/ 



(6.4) 



where denotes a 2 by 2 null matrix, and Ra are 2 by 2 diagonal matrices defined by 

Ri = diag A") 

R2 = diag A) 
R3 = diag(^(-l)"A-"'-2",(-l)"+M-"'-2n+4' 
Then by use of above F-matrices we obtain the representations of a2 as 



where is a 2 by 2 null matrix as before, and are 














Bi 







B2 




















B3 


























Bi 



















\ 








B? 





0/ 



(6.5) 



Bi 



1 



a-'^./a;;^, a-^~^4K. 



n-l . 



Bo 



A^ An + A 



n—l 



AAr. 



V An+l An-1 

\n+l A-n'^-2n+3 



V An+1 An-l 



B, 



B3 = (-l)"+M- 
We mean that B^ is a transpose of Bi. 

We note that for all a are block diagonal 

p {a,') = diag (R,2, R,2, R^2^ R^2^ R^^, R^^) 

p {a,') = diag (Bi Bj, B^ B, B^, B^ Bi, B^ B^) 
p ((T32) = diag (Ri^, R,2, R^2^ R^2^ R32, R^^) 



(6.6) 
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This fact shows that, when we restrict braiding operators generated from even powers of braid 
operators {a^,a2,cr^}, all the 2-dimensional spaces {\0a), |la)} are preserved. Then these 2- 
dimensional space may be treated as qubit, and block matrices in (|6.6I) are treated as unitary 
operators on these spaces. 

It should be noted that other sets of braid operators preserves 2-dimensional space, e.g. 

• A set of {a^, a2, cr^^} preserves \i2) and {is). 

• A set of {(Ti, (72^, 0-3} preserves l^s) and l^g). 

It is remarked that some of the braiding matrices are given in Ref. 2 based on explicit form of 
the correlation functions. 



7. Many-Quasi-Particle States and Quantum 6j Symbol 

Up to now, we have derived braiding matrices for systems of small numbers of quasi-particles 
based on the skein theory. We have derived all F-matrices in an elementary and combinatorial 
way based on the Jones-Wenzl idempotent. Although, explicit forms of the general F-matrices 
are known as the quantum 6j system, and we can construct the braiding matrices for arbitrary 
numbers of quasi-particles. 



7.1 Quantum Qj Symbol 

We borrow results on the quantum 6j symbol from Ref. 131| (see also Refs. BHS- We define 
g-integer and g-factorial by 

J^2n y4^~2n 



n 



\n-l 



nn 



ni 



The quantum 6j-symbol comes from the following Temperley-Lieb recoupling diagram; 

be b 



a b i 
c d j 



X 



(7.1) 



We have a complicated formula for the quantum 6j-symbol, which was first computed in Ref. |33 



a b i 



A, 



c d jj 9{a,d,i) 9{b,c,i) 



Tet 



a b i 
c d j 



(7.2) 



Here 9(a, b, c) is the ^^-net 



(a, b, c) 



(7.3) 
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which is given explicitly as 



with 



^(a,6,c) = (-l) 



—a + h + c 



_1 \i+j+k 



a — b + c 



a + b — c 



The tetrahedral net is defined by 



Tet 



'A 


B 


E 




C 


D 


F 


\ /d I 



which is known to be 

Tet ^ ^ ^' 
^ C D F 



[A]\[B]\[C]\[D]\[E]\[F]\ 

X 



with 



ai 
a2 



(A + D + E) 
(B + C + E) 



a^ = -{A + B + F) 
a, = ^{C + D + F) 



bi = -{B + D + E + F) 
b2 = ^{A + C + E + F) 
63 = ^ (A + B + C + D) 



As a simple application of the quantum 6j-symbol (17.11) . we have 



a a 



Note that the 6'-net appears from 



V Ac w 

^9{a,b,c) J\ 

c ' ' b b 



^ e{a,b,c) ^ 



>a4 



The quantum 6j symbol has symmetries; 

a b i 
c d j 



c d i 
a b J 
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which follows from a 180° rotation of (17.11) . Also we have 



9{a, d, i) 
A, 



a b i 
c d j 



9{c,dJ) 



j c b 
i a d 

die 
b j a 



(7.10) 



which comes from applications of the 6j-symbol to the following net; 




bed 



bed 



bed 



bed 



bed 



Figure 3: Pentagon identity 



Recursive use of the recoupling diagram (17.11) gives us the orthogonality relation; 



y-v b ^d a k _ ^ 



c d j \ \b c i 



(7.11) 



Well known is the pentagon equation, or the Biedenharn-EUiott identity, which follows from 

Fig. m-. 



E 



a i ml \b c 



b i k 



b c k] \ k c 



d e j j \d rn i j \e a m j [j a i j \d e j 
Amongst others useful identity we use [0] 



(7.12) 



^ I c d j I ^ \ d c j 



v^ia b k\ ^^,,yi (d b ^ 



(7.13) 
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which is the i?-matrix [|45|1 . To prove this, we first compute as 

b d 



b d 
a c 



a c 

On the other hand, the left hand side is transformed as 

b d 



b d 

y< 

a c 



E 

k 

E 

i,k 

E 

i,k 

E 

i,k 



a b k 

d c j 

a b k 

d c j 

a b k 

d c J 

a b k 

d c j 



E 



a b i 
c d j 



a b k 
d c j 



b d 

1\ 



d b k 



m 
m 



-1 \ d b i 



c a k 

1 \d b i 

c a k 



a c 

b d 



Combining these expressions, we obtain (|7.13l) . 



7.2 Braid Relations for 4-Quasi-Particles with Arbitrary Spins 



We generalize results in previous sections. We define the state 

a b c d 



\'^ab;cd) 



A,- 



6{a, 6, i) 6{c, d, i) 



(7.14) 



The dimension of the Hilbert space, i.e., the number of i for fixed {a, b, c, d}, is the number of z's 
such that vertices (a, b, i) and (c, (i, i) satisfy the admissible condition (|2.24l) . The orthonormality 
follows from (TT^Sl) : 

(,'^ab;cd\jab;cd) 

The braiding matrices on the space spanned by \iab;cd) can be computed by use of the quantum 
6 j -symbol (17.11) and the twist formula (12.261) . For instance, we compute as follows; 



\ { 1 ^ a b k 




k ^ 



1 jc b j 

dak 



a c b d 
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Applying (17.131) and using a normalization factor, we obtain p{cr2). In summary, we have the 
braiding matrices as 

P{cri) \iab;cd) = (Af ) ^ \iba;cd) 



P(0-2) \iab;cd) = ^ 



l A,9ia,c,j)9ib,d,j) 
Aj 6{a, b, i) 6{c, d, i) 



1 \cd fa b j 



d c i 



\jacM) (7.15) 



Pio-3) \iab;cd) = (Af ) ^ \iab;dc) 



See that the quan tum 6j-symbol with the square root term in p(o"2) coincides with the unitary 
matrix in Ref. l32i Twist operators are 



P(6'l) \iab;cd) — Ag" \iab;cd) 



and so on. 
We have 



P (fJi 0-2 0-3^ 0-2 \iab;cd) = iK"") ^ Na6;cd) 



(7.16) 



(7.17) 
(7.18) 



7.3 Braid Relations for Many-Quasi-Particles with Spin-c/2 

We study (n + 3)-quasi-particle with spin-c/2. We set 

1 



cccccccc 



\s1S2 ■■■ Sn;c) 



y 



(7.19) 



Si S2 S3 

where Sj > 0. The admissible condition (12.241) is read as 

c < Si + Sj+i < 2K — c 

Si - C< Si+i < Si + c 

Si + Sj+i + c = mod 2 
where we mean sq = Sn+i = c. Thanks to (17.81 ) the normalization factor is computed as 

IY^=l^{Si-l,Si,c) 



(7.20) 



nr=i A. 



and we have 



I I I / . I 



(7.21) 



(7.22) 



i=l 



Here the dual state (siS2 ■ ■ ■ s„; c| is defined upside down with same normalization factor. 
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Actions of the braid operators cTi can be computed in essentially same method with the previous 
sections. Results are as follows; 

p((Ti) |si---Sn;c) = (A^^)~^ |si---s„;c) (7.23) 

p(cr„+2) |Si ■ ■ ■ Sn\ C) = {XTj'^ I'Sl ■ ■ ■ Sn, c) (7.24) 

and, for 1 < z < n + 1, we have 

p[(Ti+l) \si ■■■Sn] C) 



X \si ■ ■ ■ m ■ ■ ■ Sn', c) (7.25) 

We see that 

p (cTi 0-2 ■ ■ ■ 0"n+l Crn+2 " ' " CTg (Ti) |si • • ■ s„; c) = (Aq'')"^ |si ■ ■ ■ s„; c) (7.26) 

p ((o-i CT2 ■ ■ ■ cr„+2)"'^^) |si ■ ■ ■ s„; c) = (Ao°) ^^"^^^ |si ■ ■ ■ s„; c) (7.27) 



%^ Spin- 1/2 Quasi-Particle States 

In a case of c = 1, the admissible condition (17.201) is 

l<S^ + Si+i <2K-1 

= Si — 1 or = + 1 

where we mean sq = s„+i = 1. Normalization factor is written with the ^^-net, 

9{a, b, 1) = Ai„ax(a,fe) 



(7.28) 



Note that the notation 15152-53; 1) is different from (14.121) used in 6-quasi-particle case; |010; 1), 
|012; 1), and |210; 1) respectively coincide with |00), |01), and |10) in (14.121) . Remaining two 
states, |212; 1) and |232; 1), are related to |11) and \C) in (14.121) through the F-matrix (14.131) . 

The braid matrices follows from (I7.23I) - (I7.25I) . and non-zero quantum 6j-symbols therein are 

a 1 a — 1 . 
a - 2 1 a - r " 



1 


a + ll 


1 ( a 1 a + 


1 


a + lj 


f 1 a 1 a — 


1 


a-l] 


[ A„_i A„+i 


1 


a + IJ 


1 (Aa)' 


1 


a±l] 




1 


a±lj 
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%^ Spin- 1 Quasi-Particle States 

The admissible condition (17.201) with c = 2 is read as 

2<Si + Si+i <2K -2 

(7.29) 

Si+i e {si ± 2, Si} 

and So = = 2. The normalization factor is calculated by use of the following ^^-net; 

{%2^2^ for a = 6 
A2max{a,fe) for |a - 6| = 1 
others 

States, 1^) (15.31) . and <I5.8I) . respectively correspond to |si;2), and |siS2;2). Non-zero 

6 j -symbols are; 

2a 2 2a±2l _ j 2a 2 2a + 2 
2a±4 2 2a±2| 12a + 2 2 2a 



2 a 22a + 2l d (2a 2 2a ] Asa-i A 



2a + 2 2 2a + 2j Aaa+i \2a + 2 2 2a + 2j (Asa+i)^ 

2a 2 2al r2a 2 2a±2l d 



2a ±2 2 2a f Asaii 2 2a ±2] Aaaii Asa 

2 a 2 2a -2l A2a-2 Asa+i f2a 2 2a- 2l A2a-2 A2a+2 



2a 2 2a J (A^.-if A^. 2 2a + 2j (a^j^ 

2a2 2al rf^ r2a22al A2a-2 A2a+2 - 1 



2a 2 2a-2j A2a+i I2a 2 2aj A2a-i A2a+i 

2a+2 



2a 2 2a 1 rf^ A2a-i A2a+2 (2a 2 2a + 2] 1 



2a 2 2a + 2j (A2.+i)^ A2. 2 2a J A 



i2a-l 



2a 2 2a + 2l f 2a 2 2a 
2a 2 2a-2J~\2a- 2 2 2a- 

2a 22a-2l 1 r2a 22a-2l A2a-3 A2a+i 



2a-2 2 2a- 2j A2a_i [2a-2 2 2a J (A2a-i) 



It should be remarked that, in Ref. [UJ, studied was the quantum spin Hamiltonian acting on 
state (17.191 ) with spin-1 for the SU{2)3 theory. 
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7.4 Topological Entanglement Entropy of Many-Quasi-Particle States 
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We first discuss the entanglement entropy of the 4-quasi-particle states \iab;cd) (17- 141) . We set a 
state by 

\^) = ^Pr Kl>,cd) (7.30) 
i 

where pi E C and 

E = 1 

i 

Note that the sum of i runs over such that triples {i, a, b) and {i, c, d) are admissible (|2.24l) . The 
density matrix is written as 

ah c d 



p i*>(*i Ep'pJy9(<,,6,,)fl(e,<i,,)()(a,6,j)9(c.<i,j) 

a b c d 

We first consider the entanglement between A and B when Alice owns spin-a/2 and -6/2 quasi- 
particles and Bob others; we write A = {a, b} and B = {c, d}. Alice's reduced density matrix is 
depicted as 



P{a,b} = 



PiPj 

a b 



9{a, 6, i) 9{c, d, i) 9{a, b,j) 9{c, d,j) 




d] c 



a b 

(7.32) 



9{a,b,i) ^ 

a h 

where we have used (17.81) in the second equality. Then the entanglement entropy (13.51) is computed 
from (|3.13l) as 



where we have used (17.81 ) again. The identity (17.81) indicates that the expression (17.141 ) has or- 
thonormal bases in Alice's and Bob's spaces. Then we obtain the topological entanglement en- 
tropy (13.91) as 

^5} = ElP^l'l°g^^ (7.34) 

i 

Indeed A j denotes the quantum dimension di of the quasi-particle with spin-i/2 (12.301) who inter- 
twines Alice {a, 6} and Bob {c, d}. 
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To Study a case when Alice has spin-6/2 and -c/2 quasi-particles, A = {b, c} and B 
we prepare different bases using the quantum 6j-symbol (17.11) . We have 



where 



^ab:cd 



e{k,b, c)e{k,a,d) 



a b c d 



Pk 



E 



Pi 



' Aie{b,c,k)e{a,d,k) (a b k 



Ak6{a,b,i)9{c,d,i) \c d i 



{a,d}, 
(7.35) 

(7.36) 
(7.37) 



Using 



a b i 



c d jj ]d c j 

and (|7.lb|) . we indeed see that 



b a i 



, the orthogonal identities (17.111) . and symmetries (17.91) 



J2\pk\ 



Then one sees that Alice's reduced density matrix is computed as 

P{b,c} = ^PkPj' 



Ak 



A, 



6{b, c, k) 6{a, d, k) 9{b, c,j) 6{a, d,j] 




e{b, c,k) 



b c 

y 

k 

r\ 

b c 



(7.38) 



which gives 



A, 



By the same argument, we can identify the topological entanglement entropy with 



(7.39) 



(7.40) 



which shows that we have the quantum dimension of quasi-particle which intertwines Alice and 
Bob in the expression (17.361) . 

It is straightforward to study a case when Alice has only a spin-a/2 quasi-particle; A = {a} 
and B = {b, c, d}. From (17.321) . we have Alice's reduced density matrix as 



which gives 



P{a} = 



S{a} = S{a}=^OgAa 



(7.41) 



(7.42) 
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As the entanglement entropy (13.51) satisfies (13.41) . we may expect 



SZL^ = S{aM}=iogA, (7.43) 



{a,b,c} 

We shall check this formula skein-theoretically. The reduced density matrix is given by 

ah c 



hJ 




P{a,b,c} = X] ^''^ ' ^ (7-44) 



where 



Then we have 



■' >' 6{a,b,i)9{c,d,i) 9{a,b, j) 9{c, a, j) 



Tr{a,6,c} {P{a,b,c}) = 2^ "^J1J2 " " -"^Jn-Un I [[ ^"^^ ) 

1 



The replica trick (13.81) gives (|7.43l) as we expected. 

We stress again the the topological entanglement entropy 5'^''° also depends on how to define 
A and B. 

To conclude we have shown that the topological entanglement entropy is written in terms of the 
quantum dimension. Namely in the case that the state |^) with non-Abelian quasi-particles has a 
form (|3.1I) . the entanglement entropy is modified to 

SA = J2\P^f log (^-^^ (7.46) 

which includes the topological entanglement entropy 

j 

where \pjf denotes a probability, J2j \Pjf = 1' that, in the state quasi-particle with quantum 
dimension dj plays a role of Cupid who intertwines Alice and Bob. We have obtained (|7.47l) for 
4-quasi-particle states, but a computation for many-quasi-particle states is same because recursive 



use of the identity (17.71) enables us to rewrite any braided trivalent graphs into a graph where only a 
single Wilson line intertwines Alice and Bob. Once we have such trivalent graph, the identity (17.81) 
proves that we indeed have orthogonal bases in a sense of (13.21) . So our formula (17.471) for the 
topological entanglement entropy is correct for other many-quasi-particle states. 
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8. Concluding Remarks and Discussions 



We have studied the topological properties of many-quasi-particle states in the Read-Rezayi 
state whose effective theory is the SU {2)k Chem-Simons theory. Applying the skein theory, we 
have obtained the unitary representation of the braiding matrices of quasi-particle states. These 
braiding matrices can be used as fundamental gates to construct other unitary gates. For exam- 
ple, by use of braiding matrices (14.311) of the 4-quasi-particle states in the Pfaffian state, we can 
construct the NOT gate and the Hadamard gate by [|20n 



p {9i cTi as ai) = 



Unfortunately as was shown in Refs. [iJ, uA this SU (2)2 theory is not universally computable, 
although proposed in Ref. Q to use another non-topological gate for universal computations. 
The simplest non-Abelian theory which is responsible for universal computation is the Fibonacci 
anyon model, i.e., SU{2)3 theory. The Solovay-Kitaev algorithm [1371 \4n supports that any uni- 
tary gates can be constructed from braiding matrices (|5.25l) and (15.261) of the Fibonacci anyons. 



although it was shown [|l6ll that the gate of large Fourier transformation cannot be^ealized ex- 

",HH. It 



actly. In fact several unitary gates were approximated using these braiding matrices 
may be interesting to construct fundamental gates such as NOT, Hadamard, and CNOT, using the 
braiding matrices of SU {2)k theory presented here. 

Using the braided trivalent graph as bases of Hilbert space of many-quasi-particle states, we 
have proposed a method to compute the bipartite entanglement entropy. Topological effect to the 
entanglement entropy was discussed in Refs. l35i |39l but extraction of topological part from the 
entanglement entropy is generally subtle in computations in physical models such as the quan- 
tum dimer model [jlSL |48|] . Our construction is suitable to avoid such subtlety due to that the 
skein theory or the Chem-Simons theory is purely topological gauge theory. We have indeed 
obtained exactly topological contribution to the entanglement entropy. We have shown that many- 
quasi-particle states have the topological entanglement entropy (17.471) which depends only on the 
quantum dimension of a quasi-particle intertwining Alice and Bob. 
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